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Abstract:  A  series  of  transient  non-linear  dynamic  finite-element  method  (FEM)  analyses  per¬ 
taining  to  the  interaction  of  a  single-ply  plain-woven  balanced  square  textile-fabric  armour  with 
a  spherical  steel  projectile  is  carried  out  in  order  to  compare  the  corresponding  results  obtained 
for  two  different  yarn  models:  (a)  a  solid  FEM  model  in  which  the  warp  and  weft  yarns  are  repre¬ 
sented  using  first-order  three-dimensional  solid  elements  and  (b)  a  membrane  model  in  which 
the  same  yarns  are  represented  using  second-order  membrane  elements.  The  analyses  are  car¬ 
ried  out  under  different  yarn-yarn  and  projectile-fabric  frictional  conditions  and  under  different 
far-field  boundary  conditions  applied  to  the  edges  of  the  fabric.  The  results  obtained  showed  that 
the  two  sets  of  analyses  yield  comparable  predictions  regarding  the  temporal  evolution  and  the 
spatial  distribution  of  the  deformation  and  damage  fields  within  the  fabric,  regarding  the  ability 
of  the  fabric  to  absorb  the  projectile’s  kinetic  energy  and  regarding  the  relative  contributions 
of  the  main  energy  absorbing  mechanisms.  The  work  also  confirmed  the  roles  yarn-yarn  and 
projectile-fabric  friction  play  in  the  impact  process  as  well  as  the  effect  of  the  far-field  boundary 
conditions  applied  to  the  edges  of  the  fabric. 
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1  INTRODUCTION 

Fibre-based  materials  and  material  systems  used  in 
lightweight  body  armour  applications  have  been  sub¬ 
jects  of  intense  research  and  development  efforts  over 
the  last  half-century.  These  efforts  have  further  inten¬ 
sified  since  the  terrorist  attacks  of  1 1  September  2001. 
The  most  effective  body  armour  currently  used  by  mili¬ 
tary  and  law  enforcement  personnel  is  based  on  multi¬ 
ply  fabrics  and  flexible  fibrous  composites.  These  are 
often  supplemented  by  a  hard  and  light  ceramic  strike- 
face  plate  (e.g.  5-10  mm  thick  B4C  plates)  of  which  the 
primary  role  is  to  help  blunt,  erode,  and/or  fragment 
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armour-piercing  rounds.  The  same  fibre-based  mate¬ 
rials  are  also  critical  components  in  the  more-rigid 
armour  panels  used  for  ballistic  protection  ofVIP  auto  - 
mobiles,  light  military  vehicles,  helicopters,  and  airline 
cockpit  doors. 

Over  the  last  several  decades  major  advances  have 
been  made  in  the  areas  of  development  and  fabri¬ 
cation  of  fibre-based  materials  for  ballistic  protec¬ 
tion.  Prior  to  the  early  1970s,  poly-amide  (nylon) 
fibres  dominated  these  classes  of  materials  and 
their  defining  features  were  pronounced  stress-strain 
non-linearities  and  relatively  large  strains-to-failure. 
Since  then,  new  polymeric  fibrous  materials  have 
been  developed  with  substantially  improved  strength, 
stiffness,  and  ballistic  performance.  Among  these 
high-performance  fibres  the  most  notable  are:  (a) 
poly-aramids  (e.g.  Kevlar®,  Twaron®,  andTechnora®); 
(b)  highly- oriented  poly-ethylene  (e.g.  Spectra®  and 
Dyneema®);  (c)  poly-benzobis-oxazole  (PBO)  (e.g. 
Zylon®);  and  (d)  poly-pyridobisimi-dazole  (PIPD) 
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(e.g.  M5®).  When  tested  in  tension,  all  these  materi¬ 
als  differ  significantly  from  the  nylon  fibres,  having 
very  high  absolute  stiffness,  extremely  high  specific 
strength  (a  ratio  of  strength  and  mass  density),  and 
quite  low  (<4  per  cent)  strains-to-failure.  These  fibres 
essentially  behave,  in  tension,  as  rate-independent 
linear  elastic  materials.  When  tested  in  transverse 
compression,  however,  these  fibres  are  similar  to  nylon 
and  can  undergo  large  plastic  deformation  without  a 
significant  loss  in  their  tensile  load-carrying  capacity. 
This  behaviour  is  quite  different  from  that  found  in  car¬ 
bon  or  glass  fibres,  which  tend  to  shatter  under  trans¬ 
verse  compression  loading  conditions.  Among  the 
polymeric  fibres  mentioned,  the  best  performer,  PBO, 
has  a  tensile  strength  of  ~5  GPa  (i.e.  more  than  three 
times  the  tensile -strength  of  the  strongest  ‘piano-wire’ 
steel  but  at  a  one-fifth  of  the  steel  density,  resulting  in 
an  ~  1 5- fold  higher  density-normalized  strength). 

As  mentioned,  fabrics  based  on  high-performance 
fibres  are  extensively  employed  in  a  variety  of  bal¬ 
listic  and  impact  protection  applications.  Despite 
the  fact  that  over  the  past  two  decades  there  has 
been  a  great  deal  of  work  done  on  understanding 
the  ballistic  behaviour  of  these  fabrics  using  vari¬ 
ous  analytical  and  numerical  techniques,  the  design 
of  fabric  armour  systems  remains  largely  based  on 
the  employment  of  extensive  experimental  test  pro¬ 
grams,  empiricism,  and  old  practices.  While  such 
experimental  programs  are  critical  for  ensuring  the 
utility  and  effectiveness  of  the  armour  systems,  they 
are  generally  expensive,  time-consuming,  and  involve 
destructive  testing.  Consequently,  there  is  a  contin¬ 
ual  effort  to  reduce  the  extent  of  these  experimental 
test  programs  by  complementing  them  with  the  cor¬ 
responding  computation-based  engineering  analyses 
and  simulations. 

Among  the  main  computational  engineering  anal¬ 
yses  used  to  model  ballistic  performance  of  flexible 
armour  the  following  classes  can  be  identified. 

1.  Finite- element  analyses  based  on  the  use  of 
pin-jointed  orthogonal  bars  to  represent  flexible- 
fabric  yarns.  The  most  notable  studies  falling 
into  this  category  of  analyses  are  those  per¬ 
formed  by  Roylance  and  Wang  [1],  Shim  et  al.  [2], 
Tim  et  al.  [3],  Shahkarami  et  al.  [4],  Johnson 
et  al.  [5],  and  Billon  and  Robinson  [6].  While  the 
pin-jointed  orthogonal-bars-based  finite-element 
analyses  have  proved  to  be  very  efficient  in  approxi¬ 
mating  the  dynamic  behaviour  of  woven  fabrics,  the 
discrete  nature  of  the  yarn  models  was  associated 
with  inherent  oversimplifications  that  significantly 
limited  the  predictive  capability  of  the  analyses.  In 
particular,  important  contributions  associated  with 
the  weave  architecture,  surface-finish,  and  friction 
governed  yarn-yarn  and  layer- to-layer  contacts  (in 
multi-layer  fabrics)  could  not  be  accounted  for. 


2.  More-detailed  full-blown  three-dimensional  con¬ 
tinuum  finite -element  analyses,  such  as  those  car¬ 
ried  out  by  Shockey  et  al.  [7],  Duan  et  al.  [8-11], 
and  Zhang  et  al.  [12],  etc.  While  these  analyses 
have  proved  to  be  powerful  tools  for  capturing 
and  elucidating  the  detailed  dynamic  response 
of  single-layer  fabrics,  they  are  computationally 
very  demanding  when  applied  to  practical  armour 
systems  that  typically  contain  30-50  fabric  layers 
per  ply. 

3.  Unit-cell  based  approaches  have  been  used  exten¬ 
sively  in  order  to  derive  the  equivalent  (smeared) 
continuum-level  (membrane/shell)  material  mod¬ 
els  of  textile  composites  from  the  knowledge  of  the 
meso-scale  fibre  and  yarn  properties,  fabric  archi¬ 
tecture,  and  inter-yarn  and  inter-ply  frictional  char¬ 
acteristics.  Among  the  most  notable  studies  based 
on  these  analyses  are  those  carried  out  by  Kawa- 
bata  etal.  [13-15]  who  introduced  simple  analytical 
models  to  capture  the  uniaxial,  biaxial,  and  shear 
behaviour  of  fabrics.  Furthermore,  Ivanov  and 
Tabiei  [16]  proposed  a  micro-mechanical  mate¬ 
rial  model  for  a  woven  fabric  (in  which  a  visco¬ 
elastic  constitutive  model  was  used  to  represent 
the  mechanical  behaviour  of  the  yarns)  for  use  in 
non-linear  finite-element  impact  simulations.  In 
deriving  the  material  model,  Ivanov  and  Tabiei  [16] 
considered  the  motion  of  the  yarn-crossover  point 
and  developed  a  procedure  for  determining  the 
equilibrium  position  of  this  point  under  the  applied 
unit-cell  strains.  Recently,  King  etal.  [17]  proposed 
a  new  approach  for  deriving  the  continuum-level 
material  model  for  fabrics,  based  on  the  proper¬ 
ties  of  the  yarns  and  the  weave  architecture,  which 
involves  the  use  of  an  energy  minimization  tech¬ 
nique  to  establish  the  relationship  between  the 
configuration  of  the  fabric  structure  and  the  micro¬ 
scopic  deformation  of  fabric  components.  Similar 
unit-cell  based  continuum-level  membrane/shell 
material  models  have  been  developed  by  Boisse 
etal.  [18]  and  Peng  and  Cao  [19].  Also,  Shahkarami 
and  Vaziri  [20]  proposed  a  model  similar  to  but 
simpler  than  that  introduced  by  King  et  al.  [17] 
and  provided  a  detailed  account  of  its  incorpora¬ 
tion  into  a  material-model  subroutine  that  can  be 
readily  coupled  with  commercial  dynamic- explicit 
finite- element  codes. 

4.  The  use  of  higher-order  membrane/shell  finite- 
element  analyses  to  represent  the  dynamic  res¬ 
ponse  of  the  fabric  under  ballistic  loading 
conditions  and  overcome  the  aforementioned  com¬ 
putational  cost  associated  with  the  use  of  full 
three-dimensional  finite-element  analyses  of  the 
yarn/fabric  structure.  Among  the  studies  falling 
into  this  category,  the  most  notable  is  the  one 
carried  out  by  Scott  and  Yen  [21].  While  the  use 
of  higher  order  membrane  elements  was  found 
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to  be  indeed  advantageous  computationally,  it 
was  never  fully  validated  by  comparing  its  results 
against  either  those  obtained  experimentally  or 
those  obtained  using  full  three-dimensional  finite- 
element  analyses. 

5.  As  mentioned  earlier,  full-blown  three-dimensional 
finite-element  models  used  in  the  simulations  of 
ballistic  performance  of  textile  armour  systems 
are  quite  powerful  but  become  computationally 
rather  expensive  when  applied  to  commercial 
multi-ply  flexible -armour  systems.  On  the  other 
hand,  higher-order  membrane  models  for  flexi¬ 
ble  armour  systems  are  computationally  advan¬ 
tageous  and  offer  the  possibility  of  carrying  out 
the  finite-element  analysis  of  ballistic  performance 
of  commercial  multi-ply  flexible  armour  systems. 
However,  the  validity/reliability  of  the  higher- 
order  membrane  models  was  never  fully  validated. 
Hence,  the  main  objective  of  this  work  is  to  carryout 
a  parallel  investigation  of  the  ballistic  performance 
of  a  prototypical  single-ply  plain-woven  balanced 
square  flexible-textile  armour  using  both  full¬ 
blown  three-dimensional  and  higher-order  mem¬ 
brane  models  of  the  flexible  armour. 

The  article  is  organized  as  follows.  The  details  regard¬ 
ing  the  computational  procedures  employed  in  this 
work  in  deriving  the  yarn-level  material  properties, 
fabric  architectures,  and  the  formulation  of  a  tran¬ 
sient  non-linear  dynamics  problem  pertaining  to  the 
impact  of  a  single-ply  plain  woven  balanced  fabric 
armour  with  a  spherical  projectile  are  overviewed  in 
section  2.  The  main  results  obtained  in  the  current 
work  are  presented  and  discussed  in  section  3.  The 
principal  points  and  conclusions  resulting  from  this 
work  are  listed  in  section  4. 


2  COMPUTATIONAL  PROCEDURE 

This  section  provides  a  more  detailed  description  of 
the  finite -element  analyses  that  were  used  to  study 
the  effect  of  the  impact  of  a  spherical  projectile 
onto  single-ply  plain-woven  square  flexible  textile 
armour. 

2. 1  Material  model 

In  this  paper,  the  fabric  is  modelled  at  a  yarn-level 
resolution,  i.e.  warp  and  weft  yarns  are  assumed  to 
be  made  of  continuum-type  material.  In  other  words, 
no  account  of  the  fibre  length-scale  material  effects 
within  the  yarns  is  taken.  The  yarn  material  is  con¬ 
sidered  to  be  orthotropic  (more  precisely,  transversely 
isotropic)  linear-elastic  material  up  to  the  point  of  frac¬ 
ture.  Since  the  fibres  within  yarns  are  aligned  in  the 
yarn-length  direction,  this  direction  is  associated  with 


Table  1  The  orthotropic  elastic  material  data  (GPa)  for 
yarn  [11] 
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the  highest  levels  of  stiffness.  To  account  for  the  fact 
that  yarns  are  made  of  bundled  fibres  and  that,  con¬ 
sequently,  fibres  juxtaposition  can  readily  take  place 
during  transverse  and  shear  loading,  all  shear  mod¬ 
uli  and  Poisson’s  ratios  are  set  to  very  small  values 
and,  likewise,  the  transverse  Young’s  moduli  are  taken 
to  be  a  small  fraction  of  the  longitudinal  (along  the 
yarn  length)  Young’s  modulus.  A  summary  of  the  yarn- 
level  orthotropic  linear  elastic  material  properties  for 
Kevlar- 129®  fibre-based  yarns  used  in  this  work  was 
taken  from  reference  [11]  and  is  given  in  Table  1.  One 
might  expect  that  in  the  absence  of  binding  resin,  the 
shear  moduli  of  the  yarn  should  be  lower  than  their 
values  shown  in  Table  1 .  Nevertheless,  these  values  are 
reasonable  and  they  reflect  the  effect  of  inter-fibre  fric¬ 
tion  within  the  yarns.  In  addition  to  the  linear-elastic 
orthotropic  material  properties  given  in  Table  1,  and 
in  accordance  with  the  limited- ductility  and  brittle 
nature  of  the  constituent  high-performance  Kevlar- 
129  fibres,  failure  of  the  yarns  is  assumed  to  be  strain 
controlled  (i.e.  to  take  place  when  the  strain  along  the 
yarn  length  reaches  a  critical  value,  the  failure  stain  of 
4  per  cent). 

2.2  Fabric  architecture  and  yarn  discretization 
into  finite  elements 

As  mentioned  earlier,  the  fabric  is  modelled  at  the 
yarn-level  resolution,  i.e.  the  warp  and  the  weft  yarns 
are  woven  to  produce  a  plain-woven  balanced  fabric. 
Two  types  of  finite -element  methods/discretizations 
(FEM)  are  considered:  (a)  first-order  eight-node  solid- 
brick  elements  are  used  to  represent  the  yarns 
(referred  to  the  solid  FEM  analysis)  and  (b)  second- 
order  eight-node  membrane  elements  are  used  to 
discretize  the  yarn  (referred  to  as  the  membrane  FEM 
analysis).  These  two  analyses  are  discussed  in  more 
detail  here. 

2.2. 1  Solid  FEM  model 

A  typical  three-dimensional  FEM  model  for  the  fabric 
used  in  this  work  is  displayed  in  Fig.  1  (a) .  For  both  warp 
and  weft  yarns,  a  side-truncated  sinusoidal  cross- 
sectional  area  was  assumed,  which  remains  constant 
along  the  yarn  length.  Typical  values  for  the  yarn- 
crimp  wavelength  of  1.64  mm,  the  volumetric  density 
of  the  fabric  of  600kg/m3,  and  the  fabric  thickness  of 
0.23  mm  are  assumed.  The  yarn  cross-section  is  dis¬ 
cretized  using  12  elements  and  one  yarn  wavelength 
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Fig-1  Three-dimensional  finite -element  models  of  a  plain-woven  simple-ply  balanced  textile- 
fabric  armour:  (a)  first-order  eight-node  solid  finite-element  representation  and  (b)  the 
corresponding  second-order  eight-node  membrane  representation 


is  discretized  using  ten  elements.  Typically,  a  solid 
FEM  model  contained  ~200  000  first-order  eight-node 
solid-brick  elements,  ~350  000  nodes,  and  ~1 000  000 
degrees  of  freedom. 


2.2.2  Membrane  FEM  model 

A  typical  membrane  FEM  model  used  in  this  work  is 
displayed  in  Fig.  1(b).  The  same  fabric  architecture  as 
in  the  case  of  the  solid  FEM  model  was  used.  Yarns 
were  discretized  using  second-order  eight-node  mem¬ 
brane  elements.  To  account  for  the  variable  yarn  thick¬ 
ness  along  its  cross-section,  a  ‘Nodal  Thickness’  option 
was  used  within  ABAQUS /Explicit.  Typically  the  mem¬ 
brane  FEM  model  consisted  of  ~36  000  second-order 
eight-node  membrane  elements,  ~50  000  nodes,  and 
~300  000  degrees  of  freedom. 


2.3  Impact  of  the  fabric  by  a  spherical  projectile 
problem  definition 

In  both  the  FEM  analyses  discussed,  the  same  type  of 
2. 1  g  solid  steel  spherical  projectile  with  a  4  mm-radius 
and  an  initial  velocity  of  lOOm/s  were  used.  The  pro¬ 
jectile  is  modelled  as  a  rigid  body  consisting  of  2095 
tetrahedron  elements  whose  all  degrees  of  freedom 
but  one  (the  translational  degree  of  freedom  normal  to 
the  fabric)  are  constrained.  Figure  2  shows  the  initial 
configuration  of  the  projectile-fabric  system  for  the 
solid  FEM  model.  The  corresponding  initial  configu¬ 
ration  for  the  membrane  FEM  model  is  quite  similar 
to  that  displayed  in  Fig.  2  and  is  not  shown  for  the  sake 
of  brevity. 

Yarn-yarn  and  projectile-fabric  contacts  are  mod¬ 
elled  using  a  robust  hard-contact-penalty  algori¬ 
thm  provided  by  the  ABAQUS /Explicit  finite -element 
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Fig.  2  The  initial  configurations  of  the  projectile-fabric  system  for  the  solid  FEM  analysis 


program  [22]  used  in  this  work.  Also,  a  simple  Coulomb 
friction  model  was  used  to  account  for  the  yarn- 
yarn  and  projectile-fabric  frictional  effects.  Four  fric¬ 
tional  conditions  were  considered:  (a)  both  the  yarn- 
yarn  friction  coefficient  p,y/y  and  the  projectile-fabric 
friction  coefficient  /rp/f  are  set  to  0.5;  (b)  no  fric¬ 
tion;  (c)  fiy/y  =  0.5  and  p,p/f  =  0;  and  (d)  /xy/y  =  0  and 
MP/f  =  0.5. 

Three  different  types  of  boundary  conditions  for  the 
edges  of  the  fabric  are  considered:  (a)  all  four  edges  are 
clamped;  (b)  two  opposite  edges  are  clamped  and  the 
other  two  are  free;  and  (c)  all  four  edges  are  free  except 
that  the  four  corner  nodes  are  fixed. 

To  prevent  hour-glassing  effects  that  may  arise  due 
to  the  use  of  reduced-integration  elements,  a  default 
value  of  hour-glass  stiffness  was  used.  No  instability- 
type  problems  were  encountered  and  hence  no  mass¬ 
scaling  algorithm  was  used  to  control  the  maximum 
stable  time  increment. 


2.4  Projectile-to-fabric  energy  transfer  during 
ballistic  impact 

Since  one  of  the  key  ballistic-performance  require¬ 
ments  of  an  armour  system  is  to  absorb  the  kinetic 
energy  carried  by  the  impacting  projectile,  a  brief 
overview  of  the  major  projectile-to-fabric  energy 
transfer  mechanisms  during  impact  is  presented  in 
this  section.  The  impact  of  a  projectile  with  the  fabric 
is  associated  with  the  initiation  of  several  phenomena, 
the  most  important  of  which  are:  (a)  a  resisting  force  is 
exerted  by  the  fabric  onto  the  projectile,  which  causes 
a  reduction  in  the  projectile  velocity;  (b)  at  the  same 
time,  the  fabric  is  being  deformed  and  accelerated; 


and  (c)  strain  waves  generated  in  the  impact  region 
propagate  along  the  yarns  toward  the  fabric  edges, 
where  they  are  reflected.  In  the  absence  of  any  external 
force  acting  on  the  projectile-fabric  system,  the  total 
energy  of  the  system  must  be  conserved.  In  general, 
the  energy  dissipation  due  to  projectile  deforma¬ 
tion,  fibre  intermolecular  friction,  wind  resistance,  and 
acoustic  losses  are  all  assumed  to  be  negligible.  Con¬ 
sequently,  any  loss  in  projectile  kinetic  energy  A£pk 
can  be  mainly  attributed  to  the  following  three  energy¬ 
absorbing  mechanisms:  strain  energy  acquisition  by 
the  yarns  Eys,  kinetic  energy  acquisition  by  the  yarns 
Ey k,  and  the  energy  loss  due  to  frictional  sliding  at  the 
yarn-yarn  and  projectile-fabric  contact  surfaces  E{. 
The  energy  conservation  principle  then  requires  that 

A£pk  =  Eys  +  Eyk  +  Ef  (1) 

As  clearly  shown  by  Duan  et  al.  [8-11],  the  loss  of 
projectile  kinetic  energy  A Epk  is  governed  (through 
the  three  aforementioned  energy-absorption  mech¬ 
anisms)  by  several  factors  such  as  the  material 
properties  of  the  constituent  fibres,  fabric  structure, 
boundary  conditions,  projectile  geometry,  impact 
velocity,  friction  between  the  projectile  and  the  fab¬ 
ric,  and  yarn-yarn  and  fib re-to  -fibre  friction  within  the 
fabric  itself. 

Traditionally,  the  overall  decrease  in  projectile 
kinetic  energy  (A£pk)  is  determined  using  ballistic 
impact  experiments  in  which  only  the  projectile  initial 
velocity  (uO  and  the  residual  velocity  (ur)  are  measured. 
The  overall  decrease  in  projectile  kinetic  energy  ( Ai?pk) 
is  then  defined  as 

A£pk  =  *m[of  -  ur2]  (2) 
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where  m  is  the  mass  of  the  projectile.  Recently,  Star- 
ratt  etal.  [23]  developed  a  ballistic-impact  method  for 
continuous  measurement  of  projectile  velocity  v(t). 
The  loss  of  projectile  kinetic  energy  as  a  function  of 
time  t  during  impact  is  then  defined  as 

A£pk  =  ^m[vf  -  u(tf]  (3) 

While  this  experimental  technique  may  greatly  help 
to  improve  our  understanding  of  the  ballistic  impact 
behaviour  of  fabrics,  at  the  present  time,  many  essen¬ 
tial  physical  phenomena  accompanying  projectile- 
fabric  interactions  can  only  be  revealed  using  numer¬ 
ical  analyses  and  simulations  like  those  used  in 
this  work. 

3  RESULTS  AND  DISCUSSION 

In  this  section,  selected  results  are  shown  and  dis¬ 
cussed  with  the  main  emphasis  placed  on  establishing 
if  the  computationally  more  efficient  second-order 
membrane  FEM  analyses  can  be  used  in  place  of  the 
more  accurate  yet  computationally  more  demanding 
solid  FEM  analyses  to  study  the  impact  of  flexible 
textile  armour  with  high-velocity  projectiles.  Since 
it  is  well  established  that  yarn-yarn  and  projectile- 
fabric  friction  and  the  far-held  boundary  conditions 
applied  to  the  edges  of  the  flexible  textile  armour  may 
play  important  roles  in  the  ability  of  the  armour  to 
absorb  the  kinetic  energy  of  the  projectile,  different 
friction  conditions  and  far-held  boundary  conditions 
were  considered  in  this  work.  It  should  be  noted  that 
the  work  of  Cheeseman  and  Bogetti  [24]  established 
that  not  only  the  far-held  boundary  conditions,  yarn- 
yarn  friction  and  projectile-fabric  frictions  but  also 
other  factors  such  as  material  properties  of  the  hbres 
and  yarns,  fabric  structure /architecture,  and  projectile 
geometry  and  velocity  may  affect  the  ability  of  armour 
to  defeat  the  projectile.  These  additional  factors  were 
either  not  considered  in  this  work  (e.g.  hbre/yarn 
material  properties,  fabric  structure/architecture,  and 
projectile  velocity)  or  were  analysed  but  the  results 
are  not  shown  since  they  are  consistent  with  the 
results  given  in  this  paper  (e.g.  the  effect  of  projectile, 
spherical  versus  cylindrical  geometry) . 

Since  the  key  functional  requirement  for  an  armour 
system  is  to  absorb  the  kinetic  energy  carried  by  the 
projectile,  a  quantitative  comparison  of  the  results 
pertaining  to  the  temporal  evolution  of  the  absorbed 
energies  (through  yarn  deformation  and  fracture,  fab¬ 
ric  acceleration,  and  frictional-sliding  based  energy 
dissipation)  obtained  using  the  solid  FEM  and  the 
membrane  FEM  analyses  is  carried  out  in  the  sub¬ 
sequent  sections.  The  results  obtained  for  the  case 
of  far-held  boundary  conditions  corresponding  to  all 
four  fabric  edges  being  clamped,  for  all  four  frictional 


conditions  and  for  the  spherical  shape  of  the  projectile 
is  discussed  in  the  next  section.  The  effects  of  the  type 
of  the  far-held  boundary  conditions  and  the  projectile 
(spherical  versus  cylindrical)  shape  are  then  discussed 
in  the  subsequent  two  sections,  respectively. 

3.1  Four-edges  clamped  far-held  boundary 
conditions 

The  results  presented  in  this  section  correspond  to  the 
case  of  fixed  boundary  conditions  applied  to  all  four 
edges  of  the  fabric.  As  discussed  earlier,  yarn-yarn 
friction  and  projectile-fabric  friction  are  generally 
present  in  experimental  and  held  tests  of  flexible- 
fabric  armour  and  as  shown  by  Tan  etal.  [25]  in  a  series 
of  post-impact  fabric-inspection  studies,  and  play 
an  important  role  in  absorbing  the  projectile  kinetic 
energy  (e.g.  yarn-crossover  domains  near  the  impact 
region  are  found  to  be  characterized  by  extensive  hbre 
breakage  and,  in  the  absence  of  friction,  substantial 
slippage  is  observed  between  the  warp  and  weft  yarns 
while  the  impact-induced  fabric-perforation  is  typi¬ 
cally  smaller  than  the  projectile  size).  Consequently, 
the  hrst  case  analysed  here  is  the  one  in  which  both  the 
yarn-yarn  friction  coefficient  py/y  and  the  projectile- 
fabric  friction  coefficient  pp/f  are  set  to  a  non-zero 
(=0.5)  value. 

3.1.1  Fabric  deformation  and  yarn  fracture 

Examples  of  the  temporal  evolution  of  deformation 
within  the  fabric  obtained  using  the  solid  FEM  and 
the  membrane  FEM  analyses  are  displayed  in  Figs  3 
and  4,  respectively,  where  side  and  top  views  of  the 
fabric  along  with  superimposed  contour  plots  of  the 
transverse  displacement  (the  displacement  normal  to 
the  fabric  surface)  are  shown.  A  simple  comparison 
of  the  results  displayed  in  Figs  3  and  4  reveals  that 
the  temporal  evolution  of  the  deformation  state  of  the 
fabric  is  quite  similar  in  the  two  analyses  and  can  be 
summarized  as  in  the  list. 

1.  Upto~16  (is,  the  shape  of  the  transverse-deflection 
wave  front  in  the  fabric  is  nearly  circular  and,  thus, 
essentially  identical  to  the  shape  of  the  projectile- 
fabric  contact-zone. 

2.  As  time  proceeds,  the  transverse-deflection  wave 
generated  within  the  principal  yarns  (the  yarns  that 
are  in  direct  contact  with  the  projectile)  propagates 
outward  and,  through  its  interaction  with  the  sec¬ 
ondary  yarns  (the  yarns  which  are  not  in  direct 
contact  with  the  projectile)  at  the  yarn  crossovers, 
causes  the  secondary  yarns  to  also  deflect  in  the 
transverse  direction.  Consequently,  at  ~ 30-35  (is, 
the  transverse-wave  front  begins  to  acquire  a  near 
square  shape,  with  the  square  centre  coinciding 
with  the  impact-zone  centre. 
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Fig.  3  The  temporal  evolution  of  deformation  in  the  fabric  for  the  solid  FEM  model  under  the 
yarn-yarn  and  projectile-fabric  friction-coefficient  conditions  of  0.5.  Contour  bands  corres¬ 
pond  to  different  values  of  the  transverse  displacement,  i.e.  the  displacements  normal  to 
the  fabric  surface 


3.  The  square-shaped  transverse  wave  continues  to 
propagate  towards  the  clamped  edges  of  the  fab¬ 
ric  and,  upon  reaching  the  edges  at  ~60-65  p,s,  it 
is  reflected  back  towards  the  centre  of  the  impact 
zone  and  the  wave  front  acquires  an  octagonal 
shape. 

4.  At  ~100-105qs,  the  projectile  completely  pene¬ 
trates  the  fabric  and  continues  to  move  at  a  residual 


velocity  of  26  m/s  in  the  case  of  solid  FEM  anal¬ 
ysis  and  at  a  velocity  of  28  m/s  in  the  case  of 
membrane  FEM  analysis.  In  the  case  of  the  solid 
FEM  analysis  fabric  failure  is  also  observed  at  the 
four  corners  of  the  clamped  fabric  patch.  A  similar 
fabric-corner  failure  mode  was  also  observed  in  the 
case  of  membrane  FEM  analysis  but  at  rather  (>  10 
per  cent)  higher  projectile  velocities  (the  results  not 
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Fig.  4  The  temporal  evolution  of  deformation  in  the  fabric  for  the  membrane  FEM  model  under 
the  yarn-yarn  and  projectile-fabric  friction-coefficient  conditions  of  0.5.  Contour  bands 
correspond  to  different  values  of  the  transverse  displacement,  i.e.  the  displacements  normal 
to  the  fabric  surface 


shown  for  brevity),  as  well  as  a  quite  larger  fail¬ 
ure  region  (i.e.  a  larger  number  of  failed  yarns  is 
observed  in  the  case  of  membrane  FEM  analysis 
(Fig.  4(d)). 

5.  It  should  be  noted  that  despite  the  fact  that  the 
projectile-fabric  model  has  two  vertical  planes 
of  symmetry,  the  damage  region  is  asymmetric 
in  both  types  of  analyses.  The  reason  for  this 
is  that  the  finite -element  discretization  of  the 


projectile  into  tetrahedron  elements  used  in  this 
work  does  not  possess  two  planes  of  symmetry. 
The  results  displayed  in  Figs  3  and  4  thus  also 
reveal  the  effect  of  small  geometrical  perturbations 
in  the  (spherical)  projectile  on  the  fabric  failure 
response. 

6.  Upon  its  penetration  by  the  projectile,  the  fabric 
recoils  and  the  final  shape  of  its  transversely- 
deflected  zone  becomes  near  conical. 
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The  results  displayed  in  Figs  3  and  4  and  their 
discussion  aforesaid  clearly  establish  that,  at  the 
yarn-yarn  and  projectile-yarn  friction-coefficient 
conditions  of  0.5,  the  two  FEM  analyses  yield 
comparable  results  regarding  the  temporal  evolu¬ 
tion  of  the  deformation  within  the  fabric  during 
impact.  However,  some  discrepancies  were  observed 
relative  to  the  extent  of  fabric  damage  and  the 
presence/absence  of  the  fabric-corner  failure  mode. 


Similar  observations  were  made  in  the  case  of  the 
three  other  frictional  conditions  studied  here.  For 
brevity,  only  the  results  corresponding  to  zero  yarn- 
yarn  and  projectile-fabric  friction  conditions  are 
displayed  (Figs  5  and  6)  and  discussed  here.  A 
simple  comparison  of  the  results  displayed  in  Figs 
3  and  4  with  the  corresponding  results  shown  in 
Figs  5  and  6  established  the  following  important 
observations. 


Fig.  5  The  temporal  evolution  of  deformation  in  the  fabric  for  the  solid  FEM  model  in  the  absence 
of  yarn-yarn  and  projectile-fabric  friction.  Contour  bands  correspond  to  different  values  of 
the  transverse  displacement,  i.e.  the  displacements  normal  to  the  fabric  surface 
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16ms  44ms 


Fig.  6  The  temporal  evolution  of  deformation  in  the  fabric  for  the  membrane  FEM  model  in  the 
absence  of  yarn-yarn  and  projectile-fabric  friction.  Contour  bands  correspond  to  different 
values  of  the  transverse  displacement,  i.e.  the  displacements  normal  to  the  fabric  surface 


1.  No  significant  changes  are  introduced  in  the  tem¬ 
poral  evolution  and  the  spatial  distribution  of  the 
transverse- deflection  wave  due  to  elimination  of 
the  yarn-yarn  and  the  projectile-fabric  friction. 

2.  The  structure  of  the  fabric  (including  its  deforma¬ 
tion  and  failure)  in  the  impact  region  and  in  the 
nearby  surrounding  regions  is  greatly  affected  in  the 
case  of  zero  yarn-yarn  and  projectile-fabric  friction 
conditions.  Specifically,  in  the  zero-friction  case, 
yarns  were  substantially  displaced  in  the  in-plane 


directions  away  from  the  centre  of  impact.  This 
finding  can  be  readily  explained  by  the  fact  that 
the  friction  at  the  yarn  crossovers  provides  resis¬ 
tance  to  the  relative  tangential  motion  of  the  yarns, 
while  such  resistance  is  absent  in  the  zero-friction 
case.  Consequently,  in  the  zero-friction  case  yarns 
impacted  by  the  projectile  are  pushed  outward,  a 
fewer  number  of  yarns  are  broken  and  the  pro¬ 
jectile  manages  to  penetrate  the  fabric  mainly  by 
‘wedging’  through  it. 
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The  temporal  evolution  of  deformation  and  failure 
state  of  the  fabric  for  the  remaining  two  frictional  con¬ 
ditions  were  found  to  be  similar  to  those  displayed 
in  Figs  3  to  6  and  the  structure  of  the  fabric  in  the 
impact  region  and  in  the  nearby  surrounding  regions 
was  more  affected  than  in  the  /ry/y  =  /xp/f  =  0.5  case 
and  less  affected  relative  to  the  zero-friction  case.  No 
significant  differences  in  the  fabric  structure  could  be 
established  between  the  /ry/y  =  0//u.p/f  =  0.5  case  and 
the  n,y/y  —  0.5/ p.P//f  =  0  case. 

The  effect  of  the  yarn-yarn  and  projectile-fabric 
frictional  conditions  on  the  fabric  failure  pattern  at  the 
moment  of  complete  fabric  penetration  by  the  projec¬ 
tile  is  displayed  in  Figs  7(a)  to  (d)  and  8(a)  to  (d)  for 
the  cases  of  the  solid  and  the  membrane  FEM  analy¬ 
ses,  respectively.  The  results  displayed  in  these  figures 
clearly  establish  the  following  points. 

1.  The  type  of  (solid  versus  membrane)  of  FEM  anal¬ 
ysis  used  clearly  affects  the  predicted  damage  pat¬ 
tern  of  the  fabric  for  each  of  the  four  yarn-yarn  and 
projectile-fabric  frictional  conditions.  Specifically, 
the  size  of  the  penetration  hole  and  the  number 
of  yarns  broken  under  given  frictional  conditions 


at  the  yarn-yarn  crossovers  and  the  projectile- 
fabric  contact  surface  is  larger  in  the  case  of  the 
membrane  FEM  analysis.  This  is  clearly  seen  by 
comparing  the  zero -friction  results  displayed  in 
Figs  7(b)  and  8(b)  where  it  is  seen  that  two  vertical 
principal  yarns  in  the  case  of  the  solid  FEM  anal¬ 
ysis  remain  unbroken,  Fig.  7(b),  while  all  principal 
yarns  are  broken  in  the  case  of  the  membrane  FEM 
analysis,  Fig.  8(b). 

2.  The  size  of  the  penetration  hole  and  the  number  of 
broken  yarns  are  both  increased  in  the  presence  of 
yarn-yarn  and  projectile-fabric  friction.  This  obser¬ 
vation  can  be  linked  to  the  fact  that  friction  at  yarn 
crossovers  creates  resistance  to  the  relative  motion 
between  the  contacting  principal  and  secondary 
yarns,  while  the  presence  of  projectile-fabric  fric¬ 
tion  hampers  in-plane  outward  deflection  of  the 
principal  yarns.  Consequently,  yarns  impacted  by 
the  projectile  become  loaded/strained  and  ulti¬ 
mately  break.  In  the  absence  of  friction,  yarns  are 
able  to  easily  slip  past  each  other  and  slip  pass  the 
projectile  during  the  impact.  As  a  result,  fewer  yarns 
become  strained  up  to  the  point  of  fracture.  Fur¬ 
thermore,  the  yarns  that  escaped  the  projectile  can 


J-ly/y  0.5}  ftp/f  0.0  (j-y/y  0 ■  0 ,  (j-p/f  0.5 

Fig.  7  The  effect  of  yarn-yarn  /zy/y  and  projectile-fabric  ji  p/f  friction  coefficients  on  the  failure 
pattern  in  the  fabric  for  the  case  of  the  solid  FEM  model 
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(a) 


M-y/y  0.5,  (-Xp/f  0.5 


(C) 


<d) 


J-ly/y  0.5,  Lip  J  0.0 


ii\  \  0.0,  Up  f  0.5 


(b) 


Fig.  8  The  effect  of  yarn-yarn  /xy/y  and  projectile-fabric  n-p/f  friction  coefficients  on  the  failure 
pattern  in  the  fabric  for  the  case  of  the  membrane  FEM  model 


relax  more  toward  their  initial  position,  making  the 
final  size  of  the  penetration  hole  smaller. 

3.  The  size  of  the  penetration  hole  and  the  num¬ 
ber  of  broken  yarns  are  both  reduced  in  the  case 
when  either  yarn-yarn  or  projectile-yarn  friction 
are  absent,  are  found  to  be  between  their  coun¬ 
terparts  in  the  /ry/y  =  /rp/f  =  0.5  and  /iy/y  =  /xp/f  =  0 
cases.  For  both  the  solid  and  the  membrane  FEM 
analyses,  friction  at  the  projectile-fabric  contact 
surface  appears  to  have  a  more  pronounced  effect 
than  inter-yarn  friction.  This  can  be  readily  estab¬ 
lished  by  starting  from  the  /ry/y  =  /xp/f  =  0  case, 
Figs  7(b)  and  8(b),  and  adding  projectile-fabric  and 
yarn-yarn  friction,  one  at  a  time  (Figs  7  (c)  to  (d)  and 
8(c)  to  (d),  respectively). 

3.1.2  Energy  absorption  capacity  of  the  fabric 

Temporal  evolution  of  the  three  main  components  of 
the  absorbed  energy  (i.e.  yarn  deformation  energy, 
the  fabric  kinetic  energy,  and  the  frictional  sliding 
energy),  of  a  sum  of  these  three  components,  and  of 
the  total  loss  in  the  projectile  kinetic  energy  for  the 
two  types  of  simulations,  the  yarn-yarn  and  projectile- 
fabric  frictional  coefficients  of  0.5,  are  displayed  in 


Figs  9(a)  and  10(a),  respectively.  The  corresponding 
plots  obtained  under  zero  friction  conditions  (no  fric¬ 
tional  sliding  energy  term  appears  in  this  case)  are 
displayed  in  Figs  9(b)  and  10(b).  It  should  be  noted 
that  all  the  absorbed-energy  components  displayed  in 
these  figures  are  normalized  with  respect  to  the  maxi¬ 
mum  value  of  projectile  kinetic  energy  loss  in  the  case 
at  hand.  A  simple  analysis  of  the  results  displayed  in 
Figs  9(a)  to  (b)  and  10(a)  to  (b)  reveals  the  following 
points. 

1.  No  significant  differences  (typically  <5  per  cent)  in 
the  corresponding  final  values  of  the  total  projectile 
kinetic -energy  loss  and  in  the  corresponding  final 
values  of  the  three  different  absorbed-energy  com¬ 
ponents.  More  significant  differences  are  seen  in 
temporal  evolutions  of  these  quantities.  These  dif¬ 
ferences  can  be  linked  to  the  associated  differences 
in  the  yarn-yarn  and  projectile-fabric  contact-zone 
geometries.  The  latter  differences  are  caused  by  the 
fact  that,  in  the  case  of  the  solid  FEM  analysis,  the 
side-truncated  sinusoidal  yarn  cross-section  was 
represented  explicitly  using  a  six-element  wide  fine 
mesh,  while  in  the  case  of  the  membrane  FEM 
analysis,  the  same  cross-section  was  approximated 
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Fig.  9  The  effect  of  yarn-yarn  /xy/y  and  projectile-fabric 
/tp/f  friction  coefficients  on  the  temporal  evo¬ 
lution  of  projectile’s  kinetic-energy  loss  and  on 
different  energy-absorbing  mechanisms  in  the 
fabric  obtained  using  the  solid  FEM  model: 
(a)  /xy/y  =  /tp/f  =  0.5  and  (b)  /xy/y  =  /x p/f  =  0.0 

using  a  two-element  wide  hexagon  and  a  nodal 
thickness  option.  Consequently,  the  geometry  of 
the  projectile-yarn  contact  interface  defined  in 
terms  of  slave  nodes  of  the  projectile  outer  surface 
and  the  master  element  based  surface  of  the  fabric 
is  quite  different  in  the  two  analyses. 

2.  The  sum  of  the  three  absorbed- energy  components 
is  equal  to  the  total  kinetic  energy  loss  by  the  projec¬ 
tile,  i.e.  strain  energy  of  the  fabric,  the  fabric  kinetic 
energy,  and  frictional  sliding  losses  at  yarn-yarn 
crossovers  and  projectile-fabric  contact  surfaces 
are  indeed  the  primary  modes  of  energy  absorption 
by  the  fabric. 


Fig.  10  The  effect  of  yarn-yarn  /xy/y  and  projec¬ 
tile-fabric  Hp/f  friction  coefficients  on  the  tem¬ 
poral  evolution  of  projectile’s  kinetic  energy 
loss  and  on  different  energy-absorbing  mech¬ 
anisms  in  the  fabric  obtained  using  the  mem¬ 
brane  FEM  model:  (a)  /xy/y  =  /xp/f  =  0.5  and 
(b)  /xy/y  =  /x p/f  =  0.0 

3.  Initially  (till  ~60-65  |is),  the  contributions  of  the 
yarn  strain  energy  and  the  yarn  kinetic  energy  to 
the  energy-absorption  capacity  of  the  fabric  are 
comparable. 

4.  At  longer  impact  times,  however,  the  yarn  strain 
energy  becomes  the  more  dominant  mechanism  of 
energy  absorption.  Relative  increase  in  the  strain- 
energy  contribution  continues  until  its  peak,  which 
coincides  with  the  breakage  of  yarns. 

5.  Before  yarn  breaking,  the  fabric  absorbs  between 
~80  and  90  per  cent  of  the  total  projectile  kinetic- 
energy  loss.  The  remainder  of  the  energy  absorption 
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takes  place  after  yarn  breaking,  i.e.  during  the  fabric 
perforation  process. 

6.  Past  the  point  of  yarn  breaking,  yarn  strain  energy 
undergoes  a  sharp  drop.  This  is  caused  by  its  con¬ 
version  into  fabric  kinetic  energy  and  the  frictional 
sliding  losses. 

7.  While  initial  contributions  of  the  frictional-sliding 
energy  losses  to  the  total  energy-absorption  capac¬ 
ity  of  the  fabric  are  quite  small,  during  the  fab¬ 
ric  perforation  stage,  this  mechanism  of  energy 
absorption  becomes  more  important. 

8.  The  presence  of  yarn-yarn  and  projectile-fabric 
friction  affects  the  total  energy  absorbing  potential 
of  the  fabric.  Namely,  the  projectile  kinetic-energy 
losses  corresponding  to  the  ‘Normalized  Absorbed 
Energy’  of  1.0  in  Figs  9(a)  and  10(a),  is  higher  by 
5-7  per  cent  than  their  counterparts  in  Figs  9(b) 
and  10(b). 

9.  The  presence  of  yarn-yarn  and  projectile-fabric 
friction  not  only  affects  the  frictional- sliding  energy 
loss  contribution  (the  contribution  which  is  rela¬ 
tively  small)  but  also  increases  the  contributions 
of  the  yarn  strain  energy  and  the  fabric  kinetic 
energy  to  the  fabric  energy  absorption  capac¬ 
ity.  This  finding  is  related  to  the  aforementioned 
increased  coupling  among  the  yarns  in  the  presence 
of  friction. 

To  help  further  explain  the  aforementioned  effect  of 
yarn-yarn  friction  and  projectile-fabric  friction  on  the 
energy-absorption  potential  of  the  fabric,  the  extent 
of  yarn  stressing  is  next  investigated.  Toward  that  end, 
distributions  of  the  von  Mises  equivalent  stress  in 
the  fabric  are  computed  and  displayed  in  Figs  11(a) 
and  (b)  for  the  /Xy/y  =  Pp/t  ==  0.5  and  Py  /y  =  /x  p ,y 
0  conditions,  respectively.  The  results  displayed  in 
Figs  11(a)  and  (b)  were  obtained  using  the  mem¬ 
brane  FEM  analysis.  Similar  results  (not  shown  for 
brevity)  were  obtained  using  the  solid  FEM  anal¬ 
ysis.  A  simple  analysis  of  the  results  displayed  in 
Figs  1 1  (a)  and  (b)  reveals  that  stress  is  primarily  located 
within  the  principal  yarns  and  that  a  larger  num¬ 
ber  of  yarns  are  stressed  in  the  case  when  friction  is 
present  at  the  yarn-yarn  and  projectile-fabric  contact 
surfaces. 

These  differences  are  associated  with  the  fact  that 
due  to  frictional  effects,  yarn  motion  is  more  con¬ 
strained,  leading  to  a  more  effective  load  transfer 
among  the  yarns.  This  in  turn  leads  to  more  yarns 
being  loaded  and  absorbing  energy  through  deforma¬ 
tion.  Along  the  same  yarn,  resistance  to  the  relative 
motion  of  yarns  with  respect  to  other  yarns  and  with 
respect  to  the  projectile  causes  more  yarns  to  be 
accelerated  in  the  projectile-motion  direction,  leading 
to  a  higher  kinetic-energy  contribution  to  the  fabric 
energy-absorption  capacity. 
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Fig.  11  The  effect  of  yarn-yarn  /xy/y  and  projec¬ 
tile-fabric  /Xp;,f  friction  coefficients  on  the  spa¬ 
tial  distribution  of  the  von  Mises  equivalent 
stress  within  the  fabric:  (a)  /xy/y  =  /xp/f  =  0.5  and 
(b)  /Xy/y  =  Mp/f  =  0-0 


3.1.3  Separating  the  roles  of  yarn-yarn  friction  and 
projectile-fabric  friction 

In  the  previous  section,  two  types  of  FEM  analy¬ 
ses  of  the  projectile-fabric  interactions  were  exam¬ 
ined  under  either  the  same  non-zero  yarn-yarn 
or  projectile-fabric  frictional  conditions  or  in  the 
absence  of  such  friction.  While  the  results  displayed 
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revealed  that  friction  plays  an  important  role  in  the 
deformation  and  fracture  response  of  the  fabric  and 
in  its  ability  to  absorb  the  kinetic  energy  carried  by 
the  projectile,  the  analyses  carried  out  in  the  previ¬ 
ous  section  were  not  able  to  separate  the  yarn-yarn 
and  projectile-fabric  frictional  contributions.  There¬ 
fore,  in  this  section  two  additional  frictional  condi¬ 
tions  are  considered:  (a)  zero  yarn-yarn  friction  and 
0.5  projectile-fabric  friction  coefficients,  and  (b)  zero 
projectile-fabric  and  0.5  yarn-yarn  friction  coefficient. 

Figures  12(a)  and  (b)  shows  the  temporal  evolutions 
of  the  energy-absorption  capacity  of  the  fabric  for  the 
four  frictional  conditions  considered  in  this  work  and 
for  the  two  types  of  FEM  analyses,  respectively.  In  each 
case,  the  results  are  normalized  with  respect  to  the 
overall  energy  absorbed  by  the  fabric  (i.e.  the  projectile 


Fig.  12  The  effect  of  the  four  different  yarn-yarn 
frictional  conditions  on  the  energy-absorbing 
capacity  of  the  fabric:  (a)  the  solid  FEM  model 
and  (b)  the  membrane  FEM  model 


kinetic -energy  loss)  in  the  solid  FEM  analysis  under 
yarn-yarn  /ry/y  and  projectile-fabric  n?/f  friction  coef¬ 
ficients  of  0.5.  The  results  displayed  in  Figs  12(a)  and 
(b)  clearly  established  the  following  points. 

1.  The  corresponding  predictions  regarding  the  tem¬ 
poral  evolution  of  the  absorbed  energy  are  very 
comparable  in  the  two  types  of  FEM  analyses. 

2.  Relative  to  the  zero  yarn-yarn  and  projectile- 
fabric  frictional  conditions,  the  yarn-yarn  friction 
increases  the  final  energy  absorbed  by  ~  10-12  per 
cent,  while  the  projectile-fabric  friction  increases 
the  energy  absorbed  by  ~12-15  per  cent. 

It  should  be  noted  that  while  both  yarn-yarn  and 
projectile-fabric  frictions  are  present,  the  fabric 
energy-absorption  capacity  is  increased  by  ~  18-22 
per  cent,  suggesting  that  simple  linear  superposition 
cannot  be  applied  for  the  types  of  frictions. 

The  aforementioned  effects  of  the  yarn-yarn  and 
projectile-fabric  frictions  can  be  explained  as  fol¬ 
lows. 

1.  Yarn-yarn  friction  at  the  crossovers  restricts  the 
relative  tangential  motion  of  the  yarns,  helping 
maintain  the  integrity  of  the  fabric. 

2.  Likewise,  projectile-fabric  friction  restricts  the  lat¬ 
eral  motion  of  the  principal  yarns  that  are,  con¬ 
sequently,  forced  to  remain  in  contact  with  the 
projectile. 

3.  For  both  the  cases  (1)  and  (2),  the  final  outcome 
is  that  more  yarns  are  being  loaded/strained  and 
accelerated,  which  increases  the  overall  energy- 
absorption  capacity  of  the  fabric. 

To  further  help  separate  the  effects  of  the  yarn-yarn 
and  projectile-fabric  frictions  on  the  fabric’s  energy- 
absorption  capacity,  contributions  of  the  three  energy 
absorbing  mechanisms  are  compared  for  the  (/xy/y  = 
0.5) / (/Xp/f  =  0)  case  and  the  p,y/y  =  /ip/f  =  0  case  (the 
results  not  shown  for  brevity).  It  is  found  that  that 
the  presence  of  friction  at  the  projectile-fabric  inter¬ 
faces  provides  a  substantial  contribution  to  the  fab¬ 
ric’s  energy-absorption  potential  in  the  later  stage  of 
impact.  To  rationalize  this  finding,  a  close  examina¬ 
tion  was  made  of  the  deformation  fields,  von  Mises 
stress  distributions  and  the  extent  of  fabric  damage 
(i.e.  the  number  of  yarns  broken  in  the  two  cases).  It  is 
found  that  projectile-fabric  friction  improves  the  fab¬ 
ric’s  ability  to  absorb  energy  through  several  effects: 

(a)  the  impact  load  (i.e.  the  reaction  force  measured  at 
the  reference  point  of  the  rigid  spherical  projectile)  is 
increased  and  its  peak  value  occurs  for  a  longer  time 
causing  a  larger  region  of  the  fabric  to  contribute  to 
the  energy  absorption  via  straining  and  acceleration; 

(b)  high-magnitude  von  Mises  stresses  are  present  not 
only  in  the  fabric  region  directly  underneath  the  pro¬ 
jectile  (as  is  the  case  when  projectile-fabric  friction 
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is  absent),  but  also  in  the  surrounding  regions;  and 
(c)  consequently,  a  larger  number  of  yarns  get  broken 
before  full  perforation  of  the  projectile  takes  place. 

It  should  also  be  noted  that  under  some  circum¬ 
stances,  the  largest  energy  absorption  of  the  fabric  is 
attained  when  the  projectile-fabric  friction  is  present 
but  no  friction  exist  at  the  yarn-yarn  crossover.  This 
finding  was  rationalized  as  follows:  due  to  the  absence 
of  yarn-yarn  friction,  broken  yarns  can  readily  recoil 
and  de- crimp  after  breakage,  causing  a  reduction  in 
fabric  tightness.  This  in  turn  increased  the  mobility 
of  the  unbroken  yarns  (located  below  the  projec¬ 
tile  allowing  them  to  undergo  more  deformation  and 
acquire  more  kinetic  energy  before  their  failure). 

It  should  be  recognized  that  the  values  of  the  yarn- 
yarn  and  the  projectile-fabric  frictional  coefficients 
of  0.5  are  rather  high.  Furthermore,  under  dynamic 
conditions  the  frictional  coefficients  should  be  further 
reduced.  Hence,  all  the  conclusions  made  here  regard¬ 
ing  the  effect  of  friction  pertain  to  an  ‘Upper-Bound’ 
friction  case. 


3.2  The  effect  of  far-field  boundary  conditions 

As  explained  earlier,  three  types  of  far-field  bound¬ 
ary  conditions  were  investigated:  (a)  all  four  fabric 
edges  clamped;  (b)  two  opposite  fabric  edges  clamped 
and  two  others  free;  and  (c)  the  four  fabric  corner 
nodes  fixed  and  the  remainder  of  the  fabric  edges 
free.  In  the  previous  section,  a  detailed  account  of 
the  results  is  presented  for  the  case  of  all  four  fab¬ 
ric  edges  being  fixed.  Due  to  space  limitations,  in  this 
section  only  a  brief  summary  of  the  results  obtained 
under  the  two  other  boundary  conditions  is  presented. 
The  main  finding  relative  to  the  initially  stated  objec¬ 
tive  of  this  work  is  that  for  the  other  two  boundary 
conditionsalso,  the  two  FEM  analyses  yielded  very 
comparable  results.  The  level  of  agreement  of  the 
corresponding  results  obtained  in  the  two  analyses 
were  very  similar  to  the  one  observed  in  the  previous 
section. 

In  general,  the  choice  of  far-field  boundary  condi¬ 
tions  applied  along  the  edges  of  the  fabric  is  found 
to  have  a  profound  effect  on  the  temporal  evolution 
and  spatial  distribution  of  deformation  and  damage 
within  the  fabric  and  on  the  fabric’s  energy-absorption 
capacity.  These  differences  could  be  readily  related  to 
the  fact  that  the  longitudinal  strain-stress  wave  that 
is  generated  at  the  centre  of  impact  zone  and  moves 
outward  toward  the  fabric  edges  reacts  differently  with 
the  clamped  and  free  edges.  That  is,  in  the  case  of 
clamped  edges,  the  longitudinal  waves  are  reflected 
back  toward  the  centre  of  the  impact  zone  and  the 
strain  (in  the  region  of  the  fabric  swept  by  the  reflected 
wave)  is  doubled.  In  sharp  contrast,  when  the  lon¬ 
gitudinal  wave  reaches  free  fabric  edges,  the  strain 


within  the  fabric  carried  by  the  wave  vanishes  and 
the  fabric  begins  to  move  longitudinally  toward  the 
impact  region.  Furthermore,  due  to  yarn-yarn  interac¬ 
tions,  the  free  fabric  edges  do  not  remain  straight  but 
rather,  tend  to  bow  toward  the  centre  of  the  impact 
zone.  These  differences  also  affected  propagation  of 
the  transverse  deflection  wave.  That  is,  due  to  atten¬ 
dant  larger  yarn  strains  in  the  case  of  clamped  edges, 
the  transverse  wave  tends  to  propagate  faster  along  the 
clamped  yarns  than  along  the  unclamped  yarns. 

The  type  of  boundary  conditions  applied  along  the 
fabric  edges  is  found  to  affect  significantly  the  ability 
of  the  fabric  to  absorb  the  projectile’s  kinetic  energy. 
More  specifically,  3-6  and  6-10  per  cent  more  kinetic 
energy  is  absorbed  by  the  fabric  when  all  four  edges 
are  clamped  relative  to  the  cases  when  two  edges  are 
clamped  and  no  edges  are  clamped,  respectively.  Also, 
the  role  of  friction  was  affected  by  the  type  of  bound¬ 
ary  conditions  applied.  That  is,  it  is  typically  found  that 
under  all  four  edges  free  boundary  conditions,  friction 
provides  more  contribution  to  the  energy  absorption 
through  the  frictional  sliding  mechanism.  In  sharp 
contrast,  in  the  case  of  two  edges  clamped  and  in  par¬ 
ticular,  in  the  case  of  all  edges  being  clamped,  friction 
caused  resistance  to  the  lateral  motion  of  the  yarns, 
promoting  their  deformation  and  acceleration  in  the 
direction  of  projectile  motion.  In  other  words,  fric¬ 
tion  increased  the  contributions  of  the  strain  energy 
and  fabric  kinetic  energy  to  the  energy-absorption 
potential  of  the  fabric. 


3.3  The  effect  of  projectile  geometry 

Few  selected  analyses  were  carried  out  to  establish  if 
the  type  (solid  versus  membrane)  of  FEM  analysis  sig¬ 
nificantly  affects  results  when  a  projectile  of  a  different 
geometry  is  used.  Toward  that  end,  the  4  mm  radius 
spherical  projectile  is  replaced  with  a  right  circular 
cylinder  with  equal  diameter  and  height  of  ~9  mm. 
The  two  projectiles  have  identical  volumes /weights. 
Due  to  space  limitations,  the  results  pertaining  to  the 
ballistic  impact  of  the  right  circular  cylindrical  projec¬ 
tile  with  the  fabric  is  displayed  but  is  briefly  discussed 
in  this  section. 

In  general,  the  level  of  agreement  between  the  cor¬ 
responding  results  obtained  in  the  previous  section  for 
the  spherical  projectile  is  also  seen  in  the  case  of  the 
cylindrical  projectile.  The  fabric’s  ability  to  absorb  the 
kinetic  energy  of  the  projectile  was  found  to  be  slightly 
(10-15  per  cent)  higher  in  the  case  of  the  cylindrical 
projectile.  This  observation  was  related  to  the  larger 
contribution  to  the  energy  absorption  arising  from 
friction  at  the  projectile-fabric  contact  surface.  Since 
the  projectile-fabric  contact  surface  area  is  larger  (par¬ 
ticularly  in  the  early  stage  of  impact)  in  the  case  of 
the  cylindrical  projectile,  more  yarns  are  impacted 
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and  constrained  from  moving  laterally  in  the  case 
of  the  cylindrical  projectile.  One  noticeable  differ¬ 
ence  between  the  fabric-failure  pattern  in  the  case 
of  the  spherical  and  cylindrical  projectiles  was  also 
observed.  That  is,  in  the  case  of  the  spherical  pro¬ 
jectile,  yarn  failure  was  primarily  localized  near  the 
centre  of  impact.  In  sharp  contrast,  in  the  case  of 
the  cylindrical  projectile,  failure  near  the  periphery 
of  the  projectile  was  observed.  These  differences  were 
related  to  the  observed  differences  in  the  spatial  dis¬ 
tribution  of  the  von  Mises  equivalent  stress  in  the 
two  cases.  Namely,  peaks  in  the  von  Mises  equivalent 
stress  are  found  near  the  periphery  of  the  cylindrical 
projectile. 


4  SUMMARY  AND  CONCLUSIONS 

A  set  of  two  transient  non-linear  dynamic  analyses 
pertaining  to  the  ballistic  impact  of  a  rigid  spheri¬ 
cal  projectile  with  a  plain-woven  balanced  single-ply 
square-shaped  textile-armour  is  carried  out.  The  two 
analyses  are  different  with  respect  to  the  geometri¬ 
cal  representation  of  the  yarns.  In  the  first  analysis, 
yarns  are  constructed  using  first-order  eight-node 
solid  elements,  while  in  the  second  analysis  second- 
order  eight-node  membrane  elements  are  used.  The 
objective  of  this  work  was  to  verify  if  the  computa¬ 
tionally  more  efficient  membrane  FEM  analysis  can  be 
used  in  place  of  the  more  accurately  solid  FEM  anal¬ 
ysis  and  yet  yield  reliable  results.  Based  on  the  results 
obtained  in  this  work  the  following  main  conclusions 
can  be  drawn. 

1.  The  two  sets  of  analyses  were  found  to  yield  quite 
comparable  results  under  a  variety  of  conditions 
such  as  yarn-yarn  and  projectile-fabric  frictional 
conditions  and  under  different  far-held  boundary 
conditions  applied  along  the  edges  of  the  fab¬ 
ric.  This  finding  is  particularly  true  with  respect 
to  the  temporal  evolution  and  spatial  distribu¬ 
tion  of  transverse  displacements  in  the  fabric  and 
the  overall  energy-absorbing  capacity  of  the  fabric 
armour. 

2.  In  general,  a  five-  to  six- fold  increase  in  the  compu¬ 
tational  efficiency  was  found  through  the  use  of  the 
membrane  type  FEM  analysis. 

3.  Frictional  conditions  present  at  the  yarn-yarn 
crossovers  and  at  the  projectile-fabric  contact  sur¬ 
face  significantly  affect  the  deformation  history 
of  the  fabric  and  its  ability  to  absorb  the  kinetic 
energy  of  the  projectile  and,  in  turn,  on  the  energy- 
absorption  efficiency  of  the  fabric. 

4.  Far-held  (fixed  versus  free)  boundary  conditions 
applied  along  the  edges  of  the  fabric  are  also  found 
to  have  signihcant  effect  on  the  temporal  evolu¬ 
tion  and  spatial  distribution  of  deformation  and 


damage  within  the  fabric  and,  in  turn,  on  the  ability 
of  fabric  to  absorb  the  projectile’s  kinetic  energy. 

5.  The  shape  (spherical  versus  right  circular  cylin¬ 
der)  of  the  projectile  was  not  found  to  affect  the 
computational  results  obtained  in  this  work  to  any 
signihcant  extent. 
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